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Abstract 

I We use the steepest descents method to study the integral kernel of a family of 

J> normal random matrix ensembles with eigenvalue distribution 
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£ P N (z 1 ,...,z N ) = Z- i e- N Zl 1 v«( Zi ) JJ N _,.|2 

l<i<j<N 

where V a (z) = \z\ a ,z£C and a G ]0, oo[. Asymptotic analysis with error estimates 
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1 Introduction and statement of the main result 



The investigation of non-Hermitian random matrices, whose elements are independent 
complex Gaussian variables without any constraint, began with the work of Ginibre [TO] . 
Applying the theory of Haar measure to the group GL (N, C) of N x N complex matrices, 



the joint probability distribution of the eigenvalues has shown to be given by (1.2) with 



V(z) = \z\ 2 and the eigenvalue density in the complex plane, defined by 

J pn{z)(1 2 z = — E {eigenvalues in A}) 

for any Borel set AcC, where E (•) is the expectation with respect to P/v, has shown to 
converges to the so called circular law 

pW = ( I if W - 1 ■ (1-1) 

[ otherwise 

Chau and Yue [I] have subsequently introduced ensembles of random normal matrices 
in the context of the quantum Hall problem of N electrons in a strong magnetic field, 
opening a new front of research in the area of random matrices. Since normal matrices 
are unitarily equivalent to a diagonal matrix, the probability distribution of eigenvalues 
for random normal ensembles can be achieved, exactly as in the Hermitian ensembles, by 
choosing an appropriated coordinate system that factorizes the eigenvalues contribution 
from the rest (see, respectively, Section 5.3 of [6] and [5] for the Hermitian and normal 
ensembles) . 

Normal ensembles differ from the Hermitian counterpart by the statistical dependence 
of matrix elements even for Gaussian ensembles and, most importantly, by the fact that 
their eigenvalues are generically complex. Among the usual questions concerning the 
statistics of their eigenvalues there are some related with universality that remain un- 
resolved for the normal ensembles. According to the theory of random matrices, the 
eigenvalue correlations in Hermitian, and normal ensembles as well, are given by the de- 
terminant of an integral kernel whose asymptotic behavior for large N governs their decay. 
The limit integral kernel is well known to be universal for standard models of Hermitian 
ensembles (see [6] and references therein). The scenery for normal ensembles, despite of 
certain efforts in this direction, remains undisclosed. 

The present work addresses the integral kernel of ensembles of normal matrices weighed 
by e~ NV with V depending only on the absolute value of eigenvalues. We apply the steep- 
est descents method to obtain scaling limits for the integral kernel with error estimates in 
power oil/N. Our results can be extended for a large class of radial symmetric potentials 



V satisfying condition (1.3) but we shall restrict ourselves to a sub class of potentials (1.7), 
for simplicity. Although Chau and Zaboronsky [5] have given asymptotic expressions for 
one and two-point correlation functions, the integral kernel of normal random matrices 
has not been previously considered for the models addressed here. 

The eigenvalue probability distribution of the ensemble of random normal matrices is 
given by 

P N (z ir .. ,z N ) = Z N 'e- N ^ v ^ J] \z t - Zj \ 2 (1.2) 

l<i<j<N 
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with potentials V : C — > R satisfying the properties: (i) V is continuous and (ii) 

lim (—^- - logz J = oo (1.3) 



|z|— >oo V 2 



to avoid the eigenvalues escape to infinity (see e.g. Saff and Totik [T7]). Equation (1.2) 
can be written as ^ 

Pjv(«i, • • • ,z N ) = det (K N (z h Zj))" j=1 (1.4) 
with Kn being the integral kernel 



N 

N i 



K N (z, w) = e -^ vi - z) e-^ v{w) ^ fa (z) fa (w) (1.5) 

where {fa}^ =1 is the set of the orthonormal polynomials with respect to the inner product 
(•' -)u N with weight 

dv N (z) = e- NV ^d 2 z 
and the n-point correlation function associated to P/v can be written as 

Rn ( z u ■■■ ,Zn) = det ( k n (zi, zj))l j=1 • (1.6) 

The statistics of the eigenvalues are thus given by the asymptotic behavior of the integral 
kernel. 

The main result of this paper is as follows. 
Theorem 1.1 Let 

V a (z) = \z\ a , a > (1.7) 
be a family of radially symmetric potentials, 

S(t,K) = {(eC:0<\(\<K , |argC|<^} 

be a sectorial domain of opening r and radius K and, for each < 5 < 1, let 7 = 7(0, 8) 
be such that 

cry + 8 = 1 . 



Then, the integral kernel (1.5) with V = V a satisfies 



1 ^f4^4w^"'( lw M^/« 



jy-<5+2 7 I jy 7 ' jy 7 J 47r 
with error estimation 



W(C) <0{N- & ' 2 ) (1.9) 



whenever ( G S y9/\/N, (2/a) 2 / a A^ 2 ^ 1 ^/ a \ } for some 9 > and large enough N. 
In particular, taking 8 /* 1 and, consequently, 7 \ we obtain 



1 K%{Z,W) = ^ (W) f_1 e A (ZW )* ^ ; (l + ^W)) (1.10) 



N 4vr 



TOf/i O U/Vn) error for ZW E S (e/y/N, (2/ a) 



Remark 1.2 The parameter 5 < 1 has been introduced to ensure that the eigenvalues are 
"sampling" in the bulk, out of any fixed compact domain containing the origin. The case 
of interest for applications is the limit point 8=1. The limit, as N goes to infinity, of 



any function involving the asymptotic expression (1.10) is called bulk scaling limit of 
that function. 



Remark 1.3 The restriction to a sector S(r,K) of opening r that shrinks with l/y/N 
is an artifact of our method. Equation (1.8) is expected to hold for ZW G S(t,K), 
with K = K(t) > for < r < Air /a, but our estimates on the error for replacing 
a sum by an integral, giving by the Euler-Maclaurin sum formula, break down except 
for sectors S(9N~ 13 , K) with > and (3 > 1/2 (see (4-33) and following equations). 



Numerical calculations performed in l2~Jf for a > 2 indicate that (1.8) might hold for 
ZW G S(47r/a, K) with an error decaying faster than any power of N for some K < 1 
(see also the next remark for an improved and simple estimate for a = 2). There, a 
different error: 



sup 

\z\,\w\<(2/a) 1 / a ;\a,rg(zw)\<2ir/a 



a 

47T 



(ZW) 



_1 N s 



E a /(zw) 



denoted by R", has been considered. We warn that the result of Theorem \11\ has been 
imprecisely stated in Eq. (9) of WT§ . 



Remark 1.4 Taylor remainder formula can be used to estimate the difference between 
the Taylor polynomial Sn and the function f^, respectively defined by (3.3) with 8 = 1 
(see also (4-4)) an d by the infinite sum with the same summand. For a = 2, /jv(C) = 
NQe ^ /it. By (3.2), together with the Lagrange remainder, one gets (1.10) with the error 
function satisfying |£jv(C)I = O (A r_1 / 2 (e e~ N ^~ a ^ e ^ , for some < a < 1 and large 
enough N (see calculations in Appendix We observe that (1.10\) with a = 2 holds 



with sup^g^jq |-Ejv(C)I = 0(l/vJV) for £ = ZW in a sectorial domain S(t,K) with 
K = K(r,a) > given by smallest solution f j{ e -( 1 - a ) KcoST / 2 + l = lj^] Taylor remainder 
method, together (perhaps) with some additional ingredient, may be extended for a > 2 
but it doesn't seems to work for < a < 2. 



Remark 1.5 The asymptotic behavior (1.10) for a = 2, without error estimation, was 
established in JMj. Whether n-point functions are universal for normal ensembles with 
weight e~ NV , where V(z) is a polynomial in \z\^ of positive degree with nonnegative coef- 
ficients, was addressed in |5]/. 



It follows from equations (1.10) and (1.6) that normal ensembles with the class of 



potentials V a are universal alike the Hermitian ensembles (see e.g. Subsection 5.6.1 of jH] 
and 1 13 1, for recent results): 



;i.ii) 



Corollary 1.6 Let r, Zi, . . . , z n be n + 1 complex numbers and writ 



^vK%{r,r) 



1 Wc thank an anonymous referee for suggesting the use of the Taylor remainder to estimate the error 
2. 



in (1.10 1 for a 

z Yor Hermitian ensembles, z. 



irK^ir^r) and the universal integral kernel K(z,w) in (1.12) are 



respectively replaced by xij (jrKj^(r, r)) and by the Sine function S (x — y) 



sin(x — y) 
ir(x-y) 
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Then, the following scaling limit for the n-point function 



lim 



n^oo ir n K%(r,r) n 



R%(Z 1 ,...,Z n )=det (K( Zi , Zj ))l j=1 



1.12) 



holds uniformly for r in any compact set of the open set {zGC:0< \z\ < (2/a) 1 / }, 
where 



K(z,w) 



1 (• 

— e \ 



\A \ w \' 



71 



;i.i3) 



is the integral kernel for the classical Segal-Bargmann space of entire functions. The 
bulk scaling limit (1.12) is universal in the sense that it is independent of the family of 
potentials V a . 

We shall address this and other issues related with the conformal invariance of the 
integral kernel (1.5) in a forthcoming paper [20]. Since the cancellations involved makes 
the implication of (1.12) far of being straightforward, a complete, although short, proof 
has been included in Appendix |B| 

For n = 2, (1.12) reads 



lim - 



J R n ( Z U Z 2) 



(K{ Zl , Zl )K{z 2 , z 2 ) - K{z u z 2 )K{z 2 , Zl )) 

1 



7T 



-\zi-z 2 \' 



;i.i4) 



a result already obtained for more general radial potentials (see Theorem 1 of jS]). Under 
the assumption that (1.10) can be extended to the sectorial domain SiAix/a, K) (this 
actually holds for a = 2. See Appendix [AJ , a change of variables in the integral Kernel 

by the function (Pn( Z ) = (z/\^N) 2 ^ a , which maps conformally |^ G C : \z\ < K a ^ 2 yfN^ 

into S(47r/a, K), yields 



lim ip' N (z) K% (ip N (z) , ip N (w)) p' N (w) = K(z, w) 

N— >oo 



;i.i5) 



where K.(z,w) is the integral kernel given by (1.13). This notion of universality has 
been called conformal universality in [21] . The estimates in Appendix [A] establishes the 
pointwise limit (1.15) in C x C for a = 2. 

Theorem 1.1 will be proven in Section |4j Sections [2] and [3] contain preliminary materi- 
als. The technical part of our result concerns with the error estimation of Euler-Maclaurin 
formula. Different methods needs to be employed depending on the regions considered in 
the sum. Appendix |A| estimates the Taylor remainder of (3.3) for 5 = 1 and a = 2 and 
Appendix [E] proves Corollary 1.6 



2 Ensemble of random normal matrices 

We begin with the following 

Definition 2.1 By normal ensembles we mean a probability measure 

P(M N )dM N = Z^e~ NTrV{MN) dM N (2.1) 
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on the set of N x N complex matrices Mn supported on the variety [Mn,M n ] 
invariant by unitary conjugation M N = U n M n Un ■ 



P (M N ) dM N = P[M N )dM_ 



'A' 



and 



(2.2) 



The elements = mf- + 1 < i < j < N of in the normal ensemble 

cannot be picked independently according to any product measure, absolutely continuous 
with respect to the Lebesgue measure J~| drnf^dm 1 ^ in R N +7V , even when the weight 

6 ~nttV(m n ) j g Gaussian, in view of the constraint on elements m^- with i > j ' rl So, the 
elements of Mn when sampling on normal ensembles are always statistically dependent. 
Note that the set of normal matrices with simple spectra is open and dense in IR +N and 
has full measure (see [B] for a proof in the Hermitian ensembles). 

As Mat is normal, Mn is unitarily equivalent to a diagonal matrix of eigenvalues and 
there exist Un satisfying U N X = and 



M: 



N 



(2.3) 



with Aat = diagj^i, • • • ,zn}, ordered according their absolute value: < \zj\ if % < j. 
Follow ing s ection 5.3 of [6] with few adjustments (see [5] and [8]), the spectral decompo- 
sition (2.3) considered as a change of variables M N (A^t, Un mod T N ) yields 



P(M N )dM^ 



N 



z-^e- N ^ v ^j(z, P ) n dh t n d* Pk , 

l<i<N l<j<l 



(2.4) 



where {pi} l i=1 with 21 + N = N 2 , are variables associated with the eigenvectors of M, d 2 z 
denotes the Lebesgue measure in C and 



J (z,p) 



n 

l<i<j<N 



f(p) 



is the Jacobian of ip, with / a function depending only on the eigenvectors variables 
{pi]\ = i. The eigenvalue probability distribution (2.1) of this ensemble is obtained inte- 
grating (2.4) with respect to {pi}\ =v 

The n-point correlation function is defined by (see e.g. [16] ) 



R n {zi, 



(N-n)\ 



N 



zn) 



n * 



(2.5) 



1=71+1 



and it can be written as (1.6). Stochastic processes of this form are called random deter- 
minantal point fields [18] . The present work concerns with the asymptotic analysis of the 
integral kernel (1.5) and its implications to the limit of the n-point correlation function. 
We have seen that the limit of the 2-point correlation (1.14) can be read directly from 
the asymptotic formula (1.10). The eigenvalue density p Va , associated with the normal 
ensemble defined by V a , is given by 

1 



p v - (z) 



lim i=R?(z) 



lim — K N (z, z) 



a 

47T 



z\ 



ia-2 



If they were independent, it would contradict Schur-Toeplitz statement (see e.g. 
matrix is unitarily similar to an upper (or lower) triangular matrix" . 



(2.6) 
]): "any square 
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for \z\ < {2/a) l ^ a (see Remarks 3.4, for more comment on this). Note that p Va (z)d 2 z and 
the equilibrium or extremal measure d&(z) (see e.g. [11]) agree and are supported on the 
same domain. 



3 Integral kernel of normal ensembles defined by V a 
and various estimates 



The present section is devoted to preliminary results on the integral kernel (1.5). 



Let L 2 (C, v) denote the Hilbert space of square-integrable complex-valued functions 

\f(z)\ 2 du(z) <oo 



with respect to a positive finite Borel measure v on C which, in order to ensure that all 
analytic polynomials belong to the space is assumed to satisfy 

\z\ 2n du(z) < oo , n G N . 

If P/v (C, v) denotes the iV-dimensional linear vector space of analytic polynomials of 
degree less than or equal N — 1, endowed with the inner product 



(P,0) v = / V{z)q{z)du(z) , (3.1) 
Jc 

we have 

Proposition 3.1 For each N G M, the monomials 



a 

' 1 " w 2vrr (2j/a) 



z) = J — — — —N jla z j ~ l 



with j = 1, . . . , N, form an orthonormal set in Pn (C, v%) with respect to 

dv a N (z) = e~ ma d 2 z , a > . 



The integral kernel (1.5) reads in this case 

K% (z, w) = e-%W a e-%W a K% (z, w) (3.2) 

where 

K-( ZW ) °f "'"'W" (33) 



is a reproducing kernel on Pn (C, 



Remark 3.2 For the Bergman space A 2 (Q) of square-integrable single-valued analytic 
function on a compact domain Q , there always exist a complete set of orthonormal poly- 
nomials {4>j(z)}JL 1 and the integral kernel 

oo 

K(z,w) = ^(f) j (z)(j) j (w) 

3=1 



N 



converges limjy-^oc <t>j(z)<f>j(w) = K(z,w) uniformly for any z,w in fl [2]. This is 



3=1 



not necessarily the case for unbounded domain but the same properties hold for Segal- 
Bargmann spaces A 2 (C; v) of single-valued analytic functions in C, square-integrable 
with respect to e - ' 2 ' d 2 z. We call the reader's attention to the N dependence on the inner 



product (3.1) and the fact that the limit N to infinity in (3.2) involves also a limit of the 



measure v%. As one sees from (1.10), together with 



I i a i i cc 

\z\ \w\ 



+ 



^{zw) a/2 = -\z a ' 2 -w a ' 2 \ 2 >Q 



(equality iff z = w) and equation (2.6), the limit as N — )■ oo of K% (z, w) goes for z ^ w 
and diverges for z = w. 



We shall use (3.3) to obtain an asymptotic expression as stated in Theorem 

Proof of Proposition 3.1[ We need to verify that the monomials are orthogonal with 
respect to the inner product (3.1). Writing 



4>j 0) 



with z = re ld , we have 



1* (z) , 0j (z)L» 



1 



2-Kx/h.L 



z k z J e -N\zH Z ) d 2 z 



kJ-j 



\j hlj Jo 



r k +j +l e -Nr" dr _ 



2tt 



2- 



with the Kroneker delta function S^j = 1 if k = j and otherwise, provided 

Jo a \a 

Consequently, any analytic polynomial p(z) in (C, v^) can be written as 

JV 
3=1 



with Fourier coefficients 



fa,pU = / cj )j (w)p(w)e- N ^d 2 w 



Inserting (3.5) into (3.4), gives p(z) = (^K^(z, -),pj where 

N 

Yl fa (*) fa h 



K%{z,w) 



n J^N 2 ^ a (zw 



3-1 



3=1 



(3.4) 



(3.5) 



(3.6) 



□ 

Looking for an asymptotic expansion of (3.2), a complex valued function is defined 
on the positive real line M + = (0, oo) coinciding with the summand of the integral kernel 
(3.6) on N. For fixed numbers a > 0, < 8 < 1, ( e C\ {0} and N a positive integer, let 
: K + — > C be given by 

yc w = 

and note that |<7f(a:)| = g\^\(x). 



9dx) = YWa) (3 - 7) 



Lemma 3.3 Under the above conditions on a, 5, ( and N, the real valued function 
<7|£l : K + — >• M /ias a global maximum 

g Kl (x) < maxg lc] (x) = g\ C \{x*) 
at x* = x*(a, 8, \(\ ,N)>0. For N large enough so that N > N , 

k \ 2S /a IL1 «\Ci 



with k a large universal constant, the inequality 

< x* < N 



holds and 



9\tl&) = ^Kl f ATfexp(|C|tiV 5 ) + ( 3 - 9 ) 
= SlCl f ^-T + 0f4l (3-10) 



2 ISI 4 \N 5 , 

Proof. Differentiating g\n (x) with respect to x, we have 

9' K \(x) = 9K\ (x) (log (N% |C|) - ^ (M) (3.H) 

where t/> (x) = P (or) /r (a;) is the digamma function. Since g\n (x) does not vanish and 
ijj (x) belongs to a Pick class of functions that can be analytically continued through IR + 
(see e.g. [7j), as x varies in the semi-line ip (x) increases monotonously from — oo to oo 
and the maximum of g\^\ is attained at the unique solution x = x* of 

log (*5|Cl) -§*(§*)=<>. (3.12) 
For N so large that the asymptotic expansion pQ 

oo 

Tp(y)~lo g y+--J2Bv^j (3.13) 
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of digamma function at y = N a |d can be applied (i. e., y > k where k is the constant 
mention in (3.8)), we have by (3.12) 

log (n s Id 

or equivalently, 



log -x* + + O 

a Ax 



x 



aN s \(\ 



2 * a ^ 
- = x* + - + O 
A 



which establishes (3.10). The coefficients Bij in ( |3.13 ) are the Bernoulli numbers: 

t °° ) 



n=0 



t n 

'. T 
TV. 



For ( |3.8[ ), it suffices to solve aN 5 \(\* /2 < N for N. For ( |3.9[ ), we plug ( |3.10[ ) into 
^1^1 (x*). As x* is order iV* 5 and, therefore, large enough for applying Stirling formula, 



9\C\K X 



Art id 



x* / ae 
«7r \2a;* 



ICi ATI 



1 + 



1 

AT5 



icr 



tf£ e "'KI 5 1 + 



1 

AT5 



(3.14) 



□ 



Remark 3.4 Lemma 3.3 still holds for 5 = 1 provided < |C| < (2/a) 2/a . Note that 
x* = N — a/ A + O (1/iV) < iV for \Q\ = (2/a) 2//a , which defines the domain boundary of 
the density of eigenvalues (2.6) (recall C, = ZW and \Z\ , |W| < (2/a) 1//a y ). 

The limit lim7v->.oo K^(z,w)/N calculated at zw = (/N 2 / a , given by the series 

oo 

( a /27r) ^ C i_1 /r (2j/a) , 
i=i 

converges uniformly in compact sets of C to an entire function of £ of order a/2, whose 
maximum is determined, essentially, by a single term of the series, the so called central 
index j* = J*(|CD- The next result estimates the range of indices j in (3.6) the contributes 
for its asymptotic expansion for large N. 



Lemma 3.5 Let x be a point that is at least N% logiV away from the global maximum 

(3.15) 



(3.10) of g>|£| (x), that is, 



Then 



\x - x*\ > log N. 



g K \ (x) < max (g m (x+) ,g l(l (x_)) 



where x± = x* ± N 2 log N and 



9\c\ (x±) 



jY21ogJV/( Q 2|C| 



I, n f^fN 



(3.16) 



(3.17) 
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estimates that lead to (3.14) with x± in the place of x*\ 

2x±/a 



Proof. (3.16) follows by uniqueness of the maximum value. For (3.17), we repeat the 

tl 

g\c\ ( x ±) = y 

Plugging 



«7T \ 2x± 



1 + 



1 



x ± = -\(\*N*±NnogN-- + 



into each term that appears in (3.18), yields 




2tt 



A^2 



(3.18) 



eaN s \C\ ! 



2x 



± 



a |^|«/2 N 5 / 2 2\(\ a/2 N s \N 2S 

2 log 2 log 2 N 1 1 /logiV\ 

^P | 1 T + „,i„a ^ + ^^72^5 + ( J^5j2 ) 



a|c| a/2 W a 2| Cr tf, 2|c|( 



where we have used 



= exp (1 - log(l + k)) = exp (\ - « + y + O (k 3 )^) 



and, therefore, 



><«N 5 \(\^\ 2X±/a ( ]n «/2 M s 2 2 J A ± Jlog 3 A' 

= exp |C| 7 — tttt log AT I | 1 + ' 



2x ± 



2 l/W 2 



a* |C| 



jV<V2 



Replacing in (3.18), together with (3.9), results (3.17) 



We need one more ingredient. 
Lemma 3.6 Let f : [a, b] — > M> be a convex function: 

f (Ax + (1 - X)y) < Xf ( X ) + (1-X)f(y) 
for any x, y G [a, b] and < A < land let 

P : a = x < ■ ■ ■ < Xk = b 
be the partition of [a, b] into K equally spacing subintervals of length A: 

Xj = a + jA, j E {0, . . . , K} . 
Define tj G [xj,Xj + i] by the mean value theorem: 



Xj+l 



□ 



f(x)dx = f(t j )A . 



(3.19) 
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Then, the error in the trapezoidal approximation to the integral 



s(/;^):=E / / 



(x) dx - -{f (xj) + f (x j+ i)) A 



(3.20) 



is bounded by 



0>S(/;P)> 



/(to) | | /(a*) fJtKY ]A 



(3.21) 



Proof. Without loss of generality, we suppose that / is a positive convex function. Let 



{kj} 2 j=o ^ e a numer i ca l sequence defined by 



"2j 



k 2 j+l 



Xj+l 



f (x) dx 



/ fo+i) A 



(3.22) 



for j G {0, . . . , K — 1} and note that, by the mean value theorem (3.19), 

k 2j = f(t j )A (3.23) 

for some tj G [xj, Xj + x]. We shall prove, by a geometric argument together with the 
convexity of /, that the following inequality 



h < 



h+i + k. 



i-l 



(3.24) 



holds for each i G {1, . . . , 2K — 1}. 



Since / is convex, the inequality (3.24) for i = 2j: 

-Xj+l 



f (x) dx = k 2 j < 



k 2 j+i + k 2 j-i f (x j+ i) + / (xj) 



A 



2 2 

is verified comparing the area under the function / in the interval [xj,Xj + i] (left side of 
(3.24)) with the area of a trapezoid formed by the points (xj,0) , (x J+ i,0) , (xj,f(xj)) 
and (x J+ i, / (xj+i)) (right side of (3.24)). 

Once again, by convexity of /, the inequality (3.24) for i = 2j + 1: 



f ( X 3+l) A = k 2j+l < 



k 2 j + k 2 j +2 _ 1 



f(t)dt 



(3.25) 



can be verified comparing the area under the function / in the interval [xj,x J+ 2] (2x the 
right side of (3.25 )) with the area of a rectangle of base in the interval [xj, Xj +2 ] and height 
/ (x i+i ) (2x the left side of fl3~25])). 

The later assertion is facilitate if the rectangle is replaced by a trapezoid of same area 
obtained by rotating the horizontal segment at the top around the point (xj + ±, f (x J+ i)) 
until it becomes tangent to the graph of / at that point (see figure below). 
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Now let us consider the sum 

2K 



Si = kj = k -ki-\ k 2 K-l + k 2 K 



3=0 

k i- K fU i+l - k ^ 2 -) + 'f (3.26) 

3=0 V ' 

k "- k i + £ (*« - *»- - + *«■ < 3 - 27 > 




From (3.24) and (3.26), we have 

A-l 



2" ~ ^ V i+1 2 J 



+ 



^2A . «-0 



> — + 

2 ~ 2 2 



2 A' 



and from (3.24) and (3.27), we have 

A-l 



2 



Si = fc 0-7T + ^ ( &2j 



fcjy-l + k 



2j+l 



k-iK-i , , ^ , «i k 2 K-i , , 
— 2 ^2A < k o-- 2 ™ 2K 



Since equations (3.20) and (3.27) are related by definition of {kj} 2 ^ as 



v &2A-1 . , 



the lower (3.28) and the upper (3.29) bounds yields 

h k 2K ki k 2K -i . , . , 
y + — <A: -y + S — + A; 2A - < fc 



k\ k 2K -\ , 

7, 7, r rZ 2 K 



(3.28) 



(3.29) 
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or, equivalent ly, 



kg k\ k 2 R-l 
T ~2 2 



^2K 



< E < 



which, in view of definitions (3.22) and (3.23), concludes the proof of lemma. 



□ 



Corollary 3.7 Let f : [a, b] — > K be a concave function: 

f (Ax + (1 - \)y) > Xf (x) + (1 - X)f (y) 



for all x,y G [a, b] and < A < 1 and let P, (tj). =1 and E(/;P) be as in the previous 
lemma. Then 



0< E (/;P)<|-M + M + /M_/M| 1 , 



□ 

Remark 3.8 The ideas of this proof was based in an argument used to establish the 
phenomenon of Fresnel diffraction (see e.g. JlJ/j. 



We are now in position to prove Theorem 1.1 



4 Proof of Theorem 11.1 



We shall proceed the asymptotic analysis applying the steepest descents method to the 



integral kernel (3.2). For this we assume N to be large in comparison to all other variables 



which, from now on, are kept fixed. 

It is convenient rewrite z, w and the difference of their argument using scale parameters 
7 and j3: 



Z = zN~< 
W = wN^, 7 > 



and 



9 = N 13 (arg z — arg w) , (3 > 
The equation (1.5) can thus be written as 

Z W 



N 1 W Ni 



a 
2V 



-N 1 ~ a ^\Z\ c '/2 -N 1 - a ~'\W\ a /2 



zw 



s 



N 



where 



S 



A' 



N 

E 



(]y2(l-a<y)/a Z wy 



T(2j/a) 

We introduce another auxiliary scale parameter 5 satisfying < 5 < 1 and 

cry + 5 = 1 , 



(4.1) 
(4.2) 

(4.3) 

(4.4) 

(4.5) 
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not independent. Equation (4.4) can be written as 



in order to adjust the spacing in the label that indexes the sum. Note that 7 and S are 

(4.6) 



N-l ( N 26/a ZVV y^ s 



j=0 



T (2 yj N*/a) 



where 



yj =N- 5 +jN- 5 , j = 0,...,N-l. (4.7) 

Given a function / of the class in [a, b], the Euler-Maclaurin sum formula (see 
e.g. p] with uj = and p = 1) 

AT-l 

■ 1 ' , 2 



3=0 

associated with the uniform partition P : a = y < y 1 < ■ ■ ■ < y N = b, 

Vj = a + jh, 

for j G {0, ■ • • , N — 1}, can be employed to estimate the errors 

1 



Ri = 2 (/(*) -/(«)) 



and 



i?2 



■jf (£ 1/,(G+(j+t)/i) ) rft 



in replacing the Darboux-Riemann sum of / by its integral. 
We take 



(jv M / a zif) ! 

r (2yN s /a) 



(4.9) 



l-<5 



in (4.8) with ^(x) defined by ( 3.7[ ). The partit ion N 6 = y < y\ < ■ ■ ■ < ?/jv-i = N 
of N~ s , iV 1- " 5 ] is chosen with the y/s given by (4.7). In order to simplify the notation in 
Q, from now on we fix ( = ZW = \(\ e i9/Nfi . 



Iquation (4.4) can thus be written as 

j'l V 

S N = N s 



where 



and 



/ g c (N s y) dy + ri + r 2 , 
Jn- s 



(4.10) 
(4.11) 



r 2 



-N~ 



N-l r i 

-E 

i=o 

iV-1 ^1 

-E 

3=0 



dg c (j + t) . 



(4.12) 
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The proof now proceeds in two parts. The longest one, Part I, concerns with the 
estimates of r% and r 2 . Part II applies the method of steepest descents to the integral 
term of the representation (4.8). 



I. Estimate of r x and r 2 . By the Stirling formula (see (3.14)), 



9C (AO 



(N 2S / a () 



N 




/ ae\ 2N / a 
cox \2N, 



) (iV-c) (l + 0(l/N)) = 0(N- k ) 



T (2N/a) 

holds for any power k of 1/N, in view of 2N (1 — 5) /a > 0. Since 



T (2/a) 



we conclude, by (4.11), 

r x = O (iV 25/a ) . 

According to the second mean value theorem (see e.g. [15]), for each j G {0, 
there exist G [0, 1] such that 



(4.13) 
N- 1} 



/ Q-*) d #0* + *) = -^(</c0* + 



(i)) + o(^c (i + !)-5c (i + ^)) 



Taking this into consideration, (4.12) can thus be written as 



r 2 



1 - 

2 ^ 



(4.14) 



Some considerations about (4.14) are required. We have to avoid to take absolute 
value inside the sum since any estimate that disr egards the change of sign in (4.14), leads 

r 2 to be of the leading order of the integral (4.8) given by O ( N 5 e N6 M a/2 Y\ This follows 



by (3.9) and the fact that there are 0(N S ^ 2 ) terms contributing to the sum (4.14), in view 



of Lemma 3J3 One needs to be careful and exploit the change of sign in a clever way in 
order to reduce the dependence on N from the number of terms of this sum. Because the 
estimates involve e xpon ential growth, it is convenient to divide r 2 by the maximum value 
of N s / 2 g l(l (x) (see (fi). We set 



N & / 2 g\ 



(4.15) 



for i = 1,2, and note by (4.13) that f% is exponentially smal l in N . 
Writing C = Kl e ld l Nf> with 9 G E, we have by definition (3.7) 

' ' g\ Q \ (x) cos BN~ p x + ig^\ (x) sin ON' 13 



9C [x) 



x. 



(4.16) 



As r 2 is a linear function of g^, it suffices to estimate its real part 3?e(r 2 ), since the 
estimate of Sm(r 2 ) can be done in analogous manner. 



J We have g\(\(a 



O [N 5 / 2 e N lcl ° ), the Euler-Maclaurin formula (4.101 gives an extra N s and 



N 5 / 2 results from the Gaussian integration in the steepest-descent method. See Part II. for more detail 
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The estimation of the real and imaginary parts of (4.16) depends on the period 

P= ¥\ N " 

of oscillation of {x). For this, let 11^(6) be the cardinality of the set 

A N {6) = \l G N : ^N-P < I < ^N 1 -? 

7T 7T 



(4.17) 



(4.18) 



The number Un{0) counts how many oscillations between the maximum and minimum 
value of cosON'^x there are as x varies in the interval [l,iV]. For pedagogical reason, 
we divide the estimate in two cases (i) n^(9) = 0(1) and (ii) n N {6) = 0(N £ ) for some 
< e < 1 — /3 Q The estimate for the first case can be done with less effort. In the second 
case, which may also include the previous one, the estimate is more subtle and leads to 
sharper result. 



(i) If n N (6) = n = 0(1), we write (4.14) as 



(1) , (2) 
r 2 = r 2 + r 2 



where the real part of with i = 1, 2, is given by 



Jier. 



(0 



( Ke #c u + x i)- 



(4.19) 



with A]$ being the set of points j G {1, . . . , N} such that 

> if i = 1 



^9C U + 1) - ^eg c (j + Xj 



< if i 



Let (jk)k=i denote a sequence of points right before tReg^ (j + 1) — ^Reg^ (j + Xj ) , as a 
function of j G {1, . . . , iV}, changes its sign: 



.4 
A 



(i) 

N 

(2) 
iV 



{1, . . . U {j 2 + 1, . . . , j 3 } U ■ ■ ■ U {ji-l + !,-•• , ji} 
{ji + 1, . . . , j 2 } U {j 3 + 1, . . . , j 4 } U ■ ■ ■ U {j L + 1, . . . , N} 



Since < Xj < 1 and g\(\{x) is increasing in [l,x*) and decreasing in (x*,N], the 
points (jk)t=i are essentially determined by the oscillations of the function cosON^^x 
in Jteg^ (2) = (2) cosON'^x and L = O (n^iO)) = 0(1), by hypothesis. 
By definition, we have 



(i) 



< 



(He# Ci + 1) - totfc 0')) 
|3te<? c (j! + 1) - 3f^ ? (!) + ■•• + Kep c (j L + 1) - 3fcs<fe (j L -i + 1)| 



5 We set £ = when (3 > 1. In this case rijsr(6) is always O(l). If /3 < 1, njv(#) = 0(1) when 
= 0(iV- 1 +^. 
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and 



(2) 



< 



so that 



J2 (HefocO'H-HeG/ctf + l))) 
|9f^ c {j x + 1) - Ke^ c (j 2 + 1) + • • • + Ueg c (j L + 1) - ^eg c (N + l)\ 
isp i^V^ , ^ , gK | (1) + g K | (iV + 1) 

k=l 



yields, together with ( 4.15[ ), ( |4.13[ ), Lemma 3.3 and the fact that the same holds for 
3m(r 2 ), 



(ii) Let n^{6) = 0(N £ ) for some < e < 1 — f3. Integrating (4.12) by parts gives 

N fl 



r-2 



E/ g-»)«**w+«) 



g c (j + t)dt 



£(G-)*° + < + jC 



We now split the above sum into 



r 2 = ^2 + r 2 



(4.20) 



(4.21) 



where the real part of is given by 



Ker 9 u(n) 



E (/ C? + *) - J (^C (J) + C? + !)) 

, 6A u(n) Vo ^ 



with A^-* being the set of points j G {1, • • • , iV} such that ffieg^)" (j) > (< 0). 

Let us note that the function ffieg^x) = g\Q (x) cos9x always has a well defined con- 
cavity and the cardinality of inflection points is of same order in N of the cardinality of 
critical points, since the main function responsible for both, the number of oscillations 
and changes of concavity, is the cosine. 

Let (ki) i=1 denote a sequence of points in {1, ... , N} right before (JReg^)" (j) changes 
sign. Analogously, we have 

A n 

N 



N 



{1, . . . , k x } U {k 2 + 1, . . . , k 3 } U • • • U {/c L _! + l,...,k L } 
{h + 1, . . . , k 2 } U {k 3 + 1, . . . , k 4 } U • • • U {k L + 1, . . . , N} 

where, by the same reason as in item (i), L = O (riN{6)) = 0(N e ) and, consequently, 

h+i -h = (N l - £ ) (4.22) 
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holds for i = 1, . . . , L — 1. Note also that, by ( 4.17[ ), 

e = o (jv e+/3-1 ) 



(4.23) 



Applying Lemma 3.6 (and Corollary 3.7) to each interval Jj = {ki + 1, . . . , 
z = 0, . . . , L (ko = and /c^+i = N) of size K = O (iV 1_e ) with /(a;) replaced by 3?e^(x) 
and A = 1, yields 



|9fJer 2 | < 



i=i 



Ke^ (ti) - 



Meg c (h) + Meg c (ki + 1) 



\Meg c (l)\ + \Me g( ;(N + l)\ , (4.24) 



with tk defined by the mean value theorem Meg^ (tj 



ki+l 



Meg^(x)dx. Note that the 



points (ki) i=l are closed to the inflection points (xi) i=1 of Meg^(x) and, moreover, the 
value of Meg^(x) at these points are small compared with the maximum value g^(x*). 
We shall estimate the order of Meg^(xi) and use Lemma 3.5 to reduce the number of terms 
involved in the sum (4.21). 



Taking the second derivative of the real part of (4.16), we obtain 

(Meg)" (x) = (g^ (x) - 6 2 N- 2 Pg Kl (x)) coeON^x - 29N^g\ c] (x) sinON^x 

Since derivatives of g\^\ (x) increases its value by a logarithm of N factor (see equation 
(3.11)), combined with ( |4.23[ ), it gives 



(Meg)" (x) 



(0(log 2 N) + (N 2 ^) ) cos 6N~ p x + O (N 6 ' 1 log N) sin BN~ p x . (4.25) 



0|CI W 

But we have, on the other hand 



(Meg)" (xi) = (g\' cl (x t ) - 6 2 N' 2 ^g\ C \ ( Xi )) cosfljV - ^ - 29N~ fi ^ C \ (x { ) wnON^Xi = 



holds at each inflection point x^. This together with (4.25 ) implies that the inflection point 
Xi must be at 0(1/ log N) distance from the k— th zero of cos9N~^x. Indeed, defining 
A, = 0(1/ log N) by 

(2i - 1)tt 



X j 



2\6\N-P 



+ A, 



we have 



cosON'Pxi = cos (±(i - 1/2)tt + 6N- fi A k ) = ±(-l) i sin eN~ p Ai = O (N £ ~ l / log N) 
sinON-Pxi = sin (±(i - l/2)n + ON'? Ai) =T(~iy cos 6N-?Ai = 0(1) 



and, together with (4.25), one sees that (Meg)" (x^ = holds in the leading order. Since 
the points ki, ti and ki + l are not distant from the inflection point Xi (g^ (x) varies slowly 
for each interval ki < x < ki + 1), 



|3feg c (x) 
9\c\ (x) 



\cob9N~ p x\ < O (N~ 1+£ /\ogN) 



(4.26) 



holds for x at the values {ki, U, ki + l} i=1 . 
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The number of terms that contributes to (4.21), as well as to the sum (4.24), can 
be estimated using Lemma 3.5 Instead of an interval I of size N we shall consider 



an interval I' containing x* with 0(N 5 ^ 2 log N) points. By ( 4.22 ), a number of order 
N s / 2 \ogN I N x ~ e of terms give appreciably contribution to (4.24) and, together with 



(4.26), the fact that the same estimate holds for (x) and (4.15), we conclude 

|f 2 | =0(N~ 2 ^) 
uniformly in every closed interval of < e < 1 — (5. 



II. The Method of Steepest Descents Equation (|4.10|) can be written as 

S N = N s / 2 g Kl (x* 



N S/2 



N i-s 
N~ s 



f(y)dy + h + f 2 



(4.27) 



where, by the Stirling formula (see ( |3. 14 )) 

9( {N 5 y) 



f(y) 



with 



h(y) 




N s h{y), 



0(1/N 5 )) 



ICI 



a/2 



(4.28) 
(4.29) 

Note that dteh(y) < holds for all y > and attains to its maximum §leh(y*) = at 
y* = a \C\ a ^ 2 /2 inside the domain of integration [N~ 5 , iV 1-5 ] , by condition < 5 < 1 and 
N large enough. 

We now use the steepest descents technique to estimate the integral that appears in 
(4.27). This technique uses the Cauchy theorem to deform the interval of integration 
[N~ d , N l ~ s ] into a curve C: 



2N 5 



«ICI 



a/2 



Vye NSh{y) dy 



2N & 



a |C| 



a/2 



(4.30) 



where h : C — > C is extended analytically to the complex plane, t] = y + iw and C is a 
smooth curve with extreme points rji = N~ s and 772 = iV 1-5 chosen in such a way that 
(a) it passes by the saddle point i]q = aC a ^ 2 /2 (\r)o\ = y*) defined implicitly by 



ti (Vo) 



2 , a( a/2 
a 2r] 







(4.31) 



and (b) it maximizes the function Jie/i (y, w) along a level curve 

Smli (y, w) = c 
in a neighborhood Uq of 770- If, in addition, 

9fie/i (y, w) > max{Ke/i (N~ s , 0) , M (N 1 - 5 , 0) } (4.32) 
holds along C, then the main contribution to ( 4.30[ ) will be given by the saddle point 770; if, 



on the other hand, (4.32) cannot be satisfied to any such curve C, the main contribution 



to the integral (4.30 



will be given by the extreme points. 
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At the extreme points, neither 771 nor 772 plays an important role, since both leave the 
integral (4.30) exponentially small with N. So, the contribution to (4.30) is given by the 
vicinity of the saddle point. 

Expanding h in Taylor series about r] = a( a ^ 2 /2 = a \(\ a ^ 2 e taN 06 / 2 /2, gives 



c 



a/2 



| Cr /2_ 



a 2 \(\ a/2 

with r] — 7/0 = pe %Lp G Z7o- We choose C so that 2ip — aON~^ /2 = at the the saddle point. 
Applying the steepest descents technique, the integral (4.30) can be approximate by a 
Gaussian integral in the vicinity U of r/ , resulting (see e.g. [H], for details) 

Ar ^ r /2_| f |«/2) / ^.\" / 27T?7o 



S 



N 



N s/2 g\ C \(x* 
N s/2 g\t\(x* 



1 + 



1 



^ (r / 2 _ lcr / 2) /j^r ( 1+ o(±_ 
Vicr /2 " v W* 



fi + f 2 



+ fi + f 2 



(4.33) 



Now, since by @ aiV* 5 |C| a/2 /2 < iV, 

exp (V(c /2 - icr /2 ) 



= exp |C| Q/2 (cos aN'^6/2 - 1) 
> exp (-a6 2 N 1 - 215 ) 
and, provided /3 > 1/2, it follows from the estimates of r% and r 2 in I. that 



(4.34) 



Sn = |C /2 A^exp (C^N 5 ) (1 + ^ a (C)) 



with 



|^(C)| <0(N- 5 ' 2 ) 

whenever Q G S(6N~ 1 ^ 2 , K a ' 5 ), where i'T"'' 5 = (2iV 1-5 / 'a) 2 ^ ■ Therefore, we obtain from 
(PI 



Z W 



a 

47T 



(1 + E% 6 (ZW)) (4.35) 



where we have used (4.5) with < S < 1. In particular, taking 5 /* 1, 

,,2 



a_i tv((zw) 2 -^- 



2 



(4.36) 
□ 



Remark 4.1 Equation (4-34) prevents ( = ZW to be defined in a sector S{6N~P , K a ' 6 ) 
of opening wider than 0(N~ 1 / 2 ). The introduction of the scale 5 < 1 guarantees that the 
main contribution to (4-30) comes from the saddle point for any ( 6 C fixed. Not e tha t 
K a ' s = 0(N 2 ^- S ^ a ) and for 5 = 1 we need \(\ < K a > 1 = (2/a) 2 / a (see Remark \<Q$. 
As the calculation in the appendix below indicates, \(\ may be even smaller than that, 
depending on the sector opening r. 
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A Taylor Remainder 

Let /jv(C) = N(e N ^ be a function defined for ( = \(\ e %e G C and iV a fixed natural number. 
Its Taylor remainder with respect to the polynomial Sjv(C) = + ' ' ' + jir T — r^(N() N 



of order N can be expressed by the Lagrange formula (see e.g. [19]) 



(N-l)\ 



R N (C) = f N (0 - s N (0 



9n W 



(N + l)\ 

for some < a < 1, where (?at(x) = /jvC^C); x e [0, 1] , satisfies 

g%>(x) = (rN r + N r+1 x() Ce NxC (A.l) 

for every r G N, by induction. 
Writing 

S N (() = f N (()(l + E N (0) 
the error function E N (Q = Rn(() / /jv(C) is estimated for ( in a sectorial domain S"(r, X) = 



{( G C : |arg(£)| < r/2, |£| < ^} using (A.l) together with the Stirling formula r! 
vfe(r/e) r (l + 0(l/r)): 



y/2nN 



\l + a(\ e N e-^Ci-^ICIcose^ + (l/jV)) 



so sup^ g5 ( rA -) |£?jv(C)| = O (l/^/lst^ where K = K(a,r) > is given by the smallest 
solutions of 

Ke -(l-a)K cost/2+1 = 1 ^ 

which exists and is continuous for all < a < 1 and r G [0, 2tt]. The implicit solutions of 



(A. 2) for K = K(a,r) are described in figure below for a = 1/2, 1/8 and 1/16. 

5h / 5l-i 5h 



XI— .T 

16 



2 - 




3 - 



1 - 



Oh 
2tt 



— l. oL 

In o 



2/r 
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B Proof of Corollary |1.6 



Assuming temporarily that ( |1.10[ ) holds with Z = W = r, we observe that by (1.11) 




r exp 



- + 0(l/iV) 



T + 0(1/N) 



) 



and 



arg (ZiZj) < 9/VN , 

for some 9 > and any i, j, i f N is large enough, say N > N±. We take, in addition, 
N > N where Ao is given by (3.8) with 1/ \(\ and \(\ replaced by 1/ m iny (IZiZA) and 
maxjj respectively. So, for A" > max(A* r ,A r 1 ) equation (1.10) holds with (r, r) 

and (Zj, Zj), for any in the place of (Z, W). From equation (1.9) and (1.11), it holds 



for r G C with < \r\ < (2/a) 1 / , whose closure is the support of the eigenvalues density 
(see eq. 



2.6) 



Now, applying the Taylor expansion 
(l + w) a/2 
to the exponent of (Zi,Zj), yields 



a a /a 



1 )w 2 + 0(w 3 



N[(Z l Z,) a/2 - 1 -\Z i r- 1 -\Z 3 



Aj + iVNBij + 0(1/ VN) 



(B.l) 



where 



Aij 



B 



Z{ Zj 

Xi — Xj 



2 1 Jl 



and 



A,; 



I I a/2+1 



Zj 1 

- + — ^= 

r 2\/N 



1 - 



a 



3m- 



is a real number. Let Cn and -Djv denote nx n matrices with respective entries (Cn) - = 
- exp (A i:i + iy/NBy) (l + 0(l/y/N)) and D tj = - exp {A t] ) (l + O(l/v0V)) (= with 

7T V / 7T 



-Bjj = 0). If we write A^r = diag |^exp(iv NX/) J , then CV = A N D N A N , A N A N = I (A 
and J are the complex conjugate of A^r and the identity matrix) and 

detCV = det AatDatAtv = det .DatAatAtv = detD^ ■ 



by Cauchy-Binet formula. This concludes the proof since, by (1.6) (1.10) and (B.l), the 
l.h.s of (1.12) is the determinant of a matrix whose asymptotic expansion is given by Cjv 
and 



lim det Cat = lim det = det (K (zi, zA) 



by continuity. 



□ 
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